clc;
clear;
syms x y z w v;
F=v^(1/3)*sqrt(w)*x^2*y^3*z;
tic
I=int(int(int(int(int(F,z,0,3),y,0,2),x,0,1),w,0,4),v,0,5),
toc

f=@(x)(x(1))^(1/3)*sqrt(x(2))*x(3)^2*x(4)^3*x(5);
tic,
I=quadndg(f,[0 0 0 0 0],[5,4,1,2,3]),
toc